RAL-TR-98-022 



Features of Quark and Lepton Mixing from 

Differential Geometry of Curves on Surfaces 

Jose BORDES 
bordes @ evalvx.ific.uv.es 
Dept. Fisica Teorica, Univ. de Valencia, 
■ c. Dr. Moliner 50, E-46100 Burjassot (Valencia) , Spain 

CHAN Hong-Mo 

chanhm @ v2.rl.ac.uk 

Rutherford Appleton Laboratory, 
Chilton, Didcot, Oxon, 0X11 OQX, United Kingdom 
Jakov PFAUDLER 
>• \ jakov@thphys.ox.ac.uk 

. Dept. of Physics, Theoretical Physics, University of Oxford, 

1 Keble Road, Oxford, 0X1 3NP, United Kingdom 
TSOU Sheung Tsun 
tsou @ matlis.ox.ac.uk 
00 ■ Mathematical Institute, University of Oxford, 

24-29 St. Giles', Oxford, 0X1 3LB, United Kingdom 



a: 

Oil 



Abstract 

It is noted that the CKM matrix elements for both quarks and 
(D , leptons as conceived in the Dualized Standard Model (DSM) can be 



interpreted as direction cosines obtained by moving the Darboux tri- 
hedron (a 3-frame) along a trajectory on a sphere traced out through 
changing energy scales by a 3-vector factorized from the mass matrix. 
^ \ From the 'Darboux' analogues of the well-known Serret-Frenet for- 

mulae for space curves, it is seen that the corner elements (V u b, Vtd for 
quarks, and f7 e 3, U T \ for leptons) are associated with the (geodesic) 
torsion, while the other off-diagonal elements (V us , V c d and V^, Vt s for 
quarks, and U e 2,U^\ and U^,U T 2 for leptons) with the (respectively 
geodesic and normal) curvatures of the trajectory. From this it follows 
that (i) the corner elements in both matrices are much smaller than 
the other elements, (ii) the U^, U T 2 elements for the lepton CKM ma- 
trix are much larger than their counterparts in the quark matrix. Both 
these conclusions are strongly borne out by experiment, for quarks in 
hadron decays and for leptons in neutrino oscillations, and by previous 
explicit calculations within the DSM scheme. 



The Cabibbo-Kobayashi-Moskawa matrix which gives the relative orien- 
tation in generation space between the state vectors of the three physical 
states of [/-type and D-type quarks has now been measured to a fair accu- 
racy by experiment. In the latest databook n, it is given as: 



/ 0.9745 - 0.9760 0.217- 0.224 0.0018- 0.0045 



\CKM\ 



quark 



0.217- 0.224 0.9737- 0.9753 0.036 - 0.042 
\ 0.004 - 0.013 0.035 - 0.042 0.9991 - 0.9994 

(1) 

where the modulus sign on the left-hand side means that we are giving only 
the absolute values of the entries. 

A similar matrix can in principle be defined between charged leptons 
and neutrinos, here referred to as the leptonic CKM matrix, for which the 
experimental information is still rather fragmentary in comparison, but some 
features of which are beginning to emerge. First, from the so-called muon 
anomaly observed in atmospheric neutrinos [TJ, ^ |, |], one obtains that the 
matrix element relating the muon to the heaviest neutrino z/3, namely {U^l 
is rather large, say roughly 0.45 — 0.85, according to a recent analysis [|]. 
Secondly, from the absence of any oscillation effects in reactor experiments 
such as CHOOZ [0, one obtains, at least for the squared mass difference 
Aml 3 of the two heaviest neutrinos being in the range 10~ 2 — 10~ 3 eV 2 
as favoured by the atmospheric neutrino data, that the element relating 
the electron to z^ 3 , namely \U e3 \, is small, say roughly less than 0.15 jj. 
Thirdly, for the solar neutrino puzzle ||, one is offered several solutions, 
either (i) the so-called long wave-length oscillation (LWO) solution which 
wants \U e 2\ ~ 0.4 — 0.7 |9|, |10| , or (ii) the so-called MSW solution making 



use of the Mikheyev-Smirnov-Wolfenstein [11, 12 1 mechanism which wants 



either |{7 e2 | ~ 0.4 — 0.6 (large angle solution) or |[/ e2 | ~ 0.03 — 0.05 (small 



angle solution) ||r3|j . If we choose either the LWO solution or the large angle 
MSW solution for which there seems to be some preference, we obtain for 
the leptonic CKM matrix the following tentative arrangement: 



\CKM\ 



lepton 



( * 0.4-0.7 0.0-0.15 
* * 0.45 - 0.85 I . (2) 



V 



In as much as the CKM matrix is a rotation matrix near to the identity, 
\CKM\ is roughly symmetric. And given that we are at present ignoring 
the CP-violating phase, it contains only 3 independent real parameters, so 



1 



that from the bounds of the 3 elements specified above rough bounds for the 
other elements can be calculated, but we shall not need them for our present 
discussion. 

Comparing the two matrices ([!]) and (0), one notices the following salient 
features. First, the corner elements 13 and 31 are in both cases much smaller 
than the rest. Second, the 23 element of the lepton CKM matrix is much 
larger than its counterpart in the quark CKM matrix. Third, the 12 (21) 
elements for both matrices are sizeable and of roughly the same magnitude. 
Now these features are highly significant in physical terms, for the quark 
case in explaining the branching ratios of various hadronic decays, and for 
the lepton CclS6 ; clS outlined above, in explaining certain effects in neutrino 
oscillations. It would thus be important to understand how these features 
arise in any attempt at modelling quark and lepton mixing. 

Now we have recently suggested a scheme based on the dual properties 
of the Standard Model (DSM) [jrjj which purports to explain the existence 



of exactly three generations as spontaneously broken dual colour and affords 
at the same time a means for evaluating perturbatively the CKM matrix for 
both the quark |15| and the lepton |16[ case. Calculations to 1-loop level 



have already been performed giving quite satisfactory answers which have 
the empirical features emphasized in the above paragraph. It would thus 
be interesting to find the theoretical reason within the scheme why these 
particular features should emerge. 

Let us recall some basic properties of the fermion mass matrix in the 
DSM scheme. By virtue of the manner in which the generation (dual colour) 
symmetry is broken, as suggested by the Higgs fields occurring naturally in 
the theory, the fermion mass matrix in the DSM is factorizable at tree-level, 
and even remains so after loop corrections |[L4||. Thus, 



m' = m T y' (x',y',z'). (3) 




Apart from the mass scale which may be taken to be the mass of the 
highest state of the fermion-type under consideration, the whole information 
of the mass matrix is encoded in the vector (x',y',z'). This vector rotates 
with changing energy scales, thus tracing out a trajectory in generation space, 
and in the approximation where one neglects the scale-dependence of the 
length of (x',y',z') which one cannot compute at present, one can normalize 
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the vector thus: 

x ' 2 + 2/2 + ^ = 1, (4) 

so that the trajectory lies on the unit sphere. Further, without loss of gen- 
erality we can take %' > y' > z' . 

Although the trajectory traced out by (x',y',z') can in principle be dif- 
ferent for different fermion-types it was found in ]l5j by fitting experiment 
that, to a very good accuracy, all fermion-types run on the same trajectory 
with the same speed with respect to the energy scale only differing by 
the locations of the physical states. The actual trajectory obtained in |15| 



together with the locations on it of the 12 physical states |15], [L6| are shown 
in Figure [3]. 

We shall not detail here the criterion used in the DSM for determining 
the location, namely the mass, of the lower generation states, since it is of no 
great relevance for our present discussion, but refer the reader to [jllj. What 
matters for us here is how the state vectors in generation (i.e. dual colour) 
space of the physical states are defined. 

First, for each fermion-type (i.e. whether U for U-type quarks, D for D- 
type quarks, L for charged leptons or N for neutrinos), the highest generation 
state (i.e. t,b,T,i/ 3 in the above order) is defined as that eigenstate |vi) of 
the mass matrix m! in (|3|) with non-zero eigenvalue, taken at the scale equal 
to the mass of that generation, namely at \i = mt,mb,m T ,m Ua respectively. 
It is easily seen that in Figure Q these are just the radial vectors (x' , y', z') 
themselves taken at the appropriate locations. 

Secondly, the state vector of the second generation which we called jv^) 
lSfl is defined again as the eigenvector with non-zero eigenvalue but now 



m 



of the submatrix of w! in the 2-dimensional subspace orthogonal to |vi) 
evaluated at the scale equal to the second generation mass. At the location 
on the trajectory occupied by the second generation, the vector (x',y',z') 
would have rotated to a different direction, say |vi), as illustrated in Figure 
0. The vector \v' 2 ) is the vector orthogonal to |vi) which soaks up all the 
'leakage' from this rotation. It is easy to see then that \v' 2 ) lies on the plane 
formed by |vi) and |vi) as shown in Figure |2]. Obviously, the state vector 
|vg) of the lowest generation is then, as also shown in the figure, the vector 
orthogonal to both |vi) and \v' 2 )- 

Imagine now that the locations of the highest two generations are rather 
close together, which is indeed the case for all the 4 fermion-types shown in 
Figure In that case, the vector jv^) approaches the tangent vector of the 
trajectory at (x',y',z ; ). In other words, in this approximation, the 3 state 
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Figure 1: The trajectory traced out by (x',y',z') and the locations on it of 
the 12 physical fermion states. 

vectors for the 3 generations coincide with the triad formed by the radial 
vector, the tangent vector to the trajectory, and the normal vector to both. 

Next, still in this approximation, consider the triad of state vectors for 
a pair of fermion-types, i.e. either U- with D-type quarks, or the charged 
leptons L with the neutrinos N. At the location (of the combined two high- 
est generations) of each fermion-type, we have associated a triad of state 
vectors as depicted in Figure 0. The triad at one location is obtained just 
by transporting the triad from the other location along the trajectory on the 
sphere. Now the CKM matrix is by definition the matrix giving the relative 
orientation of the triad of physical state vectors between the [/-type and D- 
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Figure 2: The state vectors of the 3 physical states belonging to the 3 gen- 
erations of each fermion-type. 



type quarks, or in the leptonic version between the charge leptons and the 
neutrinos. In view of our observations above, we conclude that it can be 
written as: 

CKM = P s expf dsA(s), (5) 

similar in form to the Dirac phase factor in gauge theory language, where 
s(/jl) is the distance on the trajectory, \i is the energy scale, and A is some 
rotation matrix yet to be specified. Notice that to 1-loop order, the distance s 
travelled per unit change in the energy scale fi is proportional to the Yukawa 
coupling strength squared p 2 , which therefore plays the role of the gauge 
coupling in the usual formula for the Dirac phase factor in gauge theory. 

To specify the matrix A, we let the two locations corresponding to the two 
fermion-types approach each other on the trajectory and evaluate the CKM 
matrix to first order in the distance As between the two locations. This can 
easily be done using some elementary results from the differential geometry 
of curves on surfaces |T7|]. Let T be a curve parametrized by arc-length s 



lying on a surface S, in our case the unit sphere. At every point on the curve 
T, we have N the unit normal to the surface, T the unit tangent to the curve, 
and we can define a third unit vector B=N A T normal to both. These form 
what is sometimes called the Darboux trihedron and are actually exactly the 
triad we are interested in, since for the sphere the radial and normal vectors 
coincide. Differentiating these vectors with respect to arc-length we get the 
following relations which we may call the Serret-Frenet-Darboux formulae 
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Figure 3: Two triads of state vectors for two fermion-types transported along 
a common trajectory. 

in analogy to the usual Serret-Frenet formulae for space curves: 



N' 

r 

B' 



KgB + K n N 



-r g N - k 9 T. 



(6) 



Here k 9 is the geodesic curvature, n n the normal curvature, and r g the 
geodesic torsion of the curve T on S. Recall that the Serret-Frenet for- 
mula for T' defines both the curvature k of the curve and its normal n at 
the given point: 

T' = ftn, (7) 

and that, being normal to T, the curvature vector kii lies on the N — B plane. 
The second equation in (^) gives its decomposition along B and N. Arranged 
in the conventional order for physicists in the form of a CKM matrix^, the 
formulae in (El) read as: 



/ 



CKM 



KgAS 



-KgAS 



\ T g As —K, n As 



-TgAS 

k„As 



(8) 



1 Notice that it is conventional to label the rows and columns of the CKM matrix 
in order of increasing mass, which is opposite to that adopted above, following |Iq ], for 
labelling the state vectors |v^) of the three generations. 
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to first order in As, this being the distance between the locations for the two 
fermion-types. This gives the matrix A in (El) as: 



/ -Kg ~Tg 



A 



K 



i) 



K n . (9) 



\ Tg K n 



Suppose now that the separation in locations between the two-fermion 
types on the trajectory is actually small as is indeed the case for the U- and 
.D-type quarks, as seen in Figure |l|. The CKM matrix is then to a good 
approximation given just by (§) above. We notice first that the off-diagonal 
corner elements of the matrix are proportional to the geodesic torsion r g while 
the other off-diagonal elements are proportional to the curvatures K g and K n . 
For the particular case of a curve on a sphere, which is our case, the geodesic 
torsion r g vanishes , so that the corner elements of the CKM matrix are of 



second order in As and therefore relatively small. Comparing this with the 
empirical CKM matrix for quarks quoted in ([[J), one sees that this prediction 
is very clearly borne out. Secondly, the elements V c b and V ts are proportional 
to the normal curvature K n which for the sphere is constant, independent 
either of the trajectory or of the location. It follows then from (^j) that their 
values are approximately given just by the distance, measured in units of 
the radius of the sphere, between the locations of t and b on the trajectory. 
From Figure ^ one estimates this to be of order of a few percent. Again, 
the prediction is seen to be borne out in ([I]). Both these special features are 
also what one finds in the explicit calculation carried out in |jl5| . As for the 
remaining off-diagonal elements V us and V c d, namely the Cabibbo angle, it 
is in a sense a little special. Being proportional to the geodesic curvature K g 
which depends both on the location and on the trajectory, it can be sizeable 
even for small separation. This is seen to be indeed the case in the empirical 
matrix (P and agrees with the fact that in the calculation [15] , the Cabibbo 



angle can be fitted to the experimental value by choosing the trajectory. 

What would happen in the case of the leptons where the separation be- 
tween the two fermion-types, i.e. between r and z^, is not so small, as seen 
in Figure [I]? We shall need in that case to evaluate the 'Dirac phase factor' 
(||) above, as was done numerically in [16|. However, even without an ex- 



plicit calculation we can already draw some qualitative conclusions. First, 
the corner element U e ^ being of second order in the distance travelled will 
remain small. This checks with the empirical estimates quoted in and 
with the calculation ||16f. Secondly, the normal curvature K n being constant 
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all along the trajectory, its contribution to the element will accumulate 
over the larger distance here and become much larger than its counterpart 
V c b in the quark case. In fact, if one just crudely takes the formula (|S|) and 
assumes that the element is still roughly proportional to the distance trav- 
elled from r to u 3 , one obtains from Figure Q that would be some 20 
times larger than V c b in the quark matrix, which estimate is seen already to 
be in fair agreement both with the empirical matrices in (H) and (|2|), and 
with the explicit calculations in |15|, |16| . Thirdly, in contrast, the remaining 
off-diagonal element U e 2 depends on the geodesic curvature n g which varies 
with the location. Its value therefore cannot be estimated from just the sep- 
aration between the locations of the two fermion-types. There is thus no 
reason to expect a value for it very different in magnitude from the values of 
V us or V c d in the quark case, nor are the values as seen in (|l|) and (Q) very 
different empirically. 

There is a further qualitative conclusion that one can draw from the 
preceding analysis about the CKM matrix which, however, cannot be checked 
against experiment but only against the calculation in fL5| . As mentioned 



above, the normal curvature n n on a sphere does not depend on the trajectory. 
Hence, the elements V c b and Vt s , being to first order in As proportional to the 
normal curvature K n , is also expected to be approximately independent of 
the trajectory. We then see why we were able, when fitting the CKM matrix 
15], first to fit the elements V c b and V ts by adjusting the Higgs vev y, and 



m 



then to fit the rest by adjusting the other vev z, for the first two elements 
will remain roughly constant in changing z which has the effect mainly of 
changing the trajectory. 

The above analysis was done in the approximations (a) that the locations 
on the trajectory of the two highest states of each fermion-type are close 
together, and (b) that the vector (x', y', z') has a normalization independent 
of the energy scale \i and hence moves on a sphere. However, the qualitative 
features deduced above for the quark and lepton CKM matrices will still re- 
main even if one relaxes these approximations so long as the scale-dependence 
of the normalization of (x',y',z') is not too strong. In any case, the CKM 
matrix can always be regarded as a 'Dirac phase factor' of the form ([5]), only 
in general the matrix A appearing there will be more involved. 

Apart from the elegance in giving us a geometrical significance for the 
CKM matrix, and a qualitative understanding of the values of its elements 
in simple geometrical terms, the above observations have the practical value, 
we think, in giving us also a nontrivial check on the basic structure of the 
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DSM scheme. The fact that the three off-diagonal elements 13, 23, and 12 are 
predicted to have quite different salient features all agreeing with experiment, 
depends crucially on the assignment in the scheme of the highest, middle and 
lowest generation states to respectively the radial vector, the tangent vector 
and the vector normal to both, on the trajectory And this trajectory is traced 
out by the vector (x f , y', z'), which is in turn obtained from a factorizable mass 
matrix. Quite a number of structural details of the scheme are thus involved 
in yielding these predictions. Further, the fact that rough values of the 
empirical matrix elements can already be deduced just from the Darboux- 
Serret-Frenet formulae @ makes the agreement with experiment obtained 
in 



lq| and [16| much less likely to be mere accidents of the calculations. 
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